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Abstract 

We construct a slowly varying space-time dependent holographic superfluid and compute 
its transport coefficients. Our solution is presented as a series expansion in inverse powers 
of the charge of the order parameter. We find that the shear viscosity associated with the 
motion of the condensate vanishes. The diffusion coefficient of the superfluid is continuous 
across the phase transition while its third bulk viscosity is found to diverge at the criti- 
cal temperature. As was previously shown, the ratio of the shear viscosity of the normal 
component to the entropy density is l/47r. As a consequence of our analysis we obtain an 
analytic expression for the backreacted metric near the phase transition for a particular type 
of holographic superfluid. 
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1 Introduction and summary 



It is quite remarkable that black holes in asymptotically AdS spacetimes can be applied, 
through the AdS/CFT correspondence [H [21 13], to myriad physical systems. Among these 
systems are gauge-theory plasmas whose hydrodynamic behavior is associated with the to- 
pography of the black hole horizon [1], and superfluid states of large N gauge-theories [SIEIE]- 
In this work we show how, in the presence of a charged scalar condensate, the black hole 
horizon encodes superfluid properties of the dual gauge theory and compute various generic 
properties of the superfluid. 

The membrane paradigm [5] suggests a relationship between fluid dynamics and black 
holes. In the context of AdS/CFT, ref. [9] made this relation precise by computing the ratio 
of the shear viscosity, r], to entropy density, s, of A/" = 4 Super- Yang Mills 

--^ (1) 

s Att 

using a Kubo formula. The ratio in ([1]) is considered to be universal in the sense that 
it is valid, in the supergravity limit, for a large class of theories [IDl CH [El [13] E The 
work of |1] provided a new layer of development in describing fluid flow using AdS/CFT by 
demonstrating a one-to-one correspondence between solutions of relativistic Navier-Stokes 
equations in field theory and solutions of Einstein's equations in the dual black hole geometry. 
Indeed, using the formalism of |1] new transport coefficients were unveiled [IHl [HI [ITJ [181 
[ini [20] and relations similar to ([T]) were found for higher order transport coefficients [20] . 

Superfluidity, initially discovered in liquid helium, can be thought of as a fluid with a 
spontaneously broken global symmetry. At temperatures below Tq ~ 2.17° K liquid helium 
undergoes a phase transition into a superfluid state. Many of the remarkable properties 
of superfluid helium can be attributed to the lack of viscosity of its superfluid component. 
In [21] and later in [22] Tisza and Landau formulated a hydrodynamic model in which a 
superfluid is described by a two component fluid: a condensate which, roughly speaking, has 
no viscosity and a normal component which is viscous. 

As we review in section [21 one can extend the Tisza-Landau model to relativistic super- 
fluids. Several relativistic versions of the inviscid superfluid can be found in the literature 
[231 [211 [23 [211 [23 [21] • viscous corrections to these models have been treated, to a certain 
extent, in [2S1 [301 [SB [32] and applied to studies of neutron stars and cold relativistic super- 



See |14j for a recent example where the authors argue that dU is, in some sense, violated in a non-isotropic 
holographic superfluid. 
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fluids in [321 133] . The number of transport coefficients associated with viscous corrections 
can be quite large. For example, the author of [M] counts 13 transport coefficients in the 
non-relativistic theory when the relative superffuid velocity is not small. 

One can make several assumptions that reduce the number of transport coefficients to 
three. Because superffuidity is typically lost if the relative velocity between the superffuid and 
normal components becomes too great, a common assumption is that this relative velocity 
is small, n'^ ^ 1. Another common assumption is the absence of parity breaking terms in 
the hydrodynamic expansion. In view of our applications to holographic superffuids with 
a traceless stress-energy tensor, we also assume conformal invariance. At leading order in 
the hydrodynamic gradient expansion, these assumptions reduce the number of transport 
coefficients to four: the shear viscosity t], a shear viscosity associated with the superffuid 
component rjs, a diffusivity k, associated with charge transport or heat conductivity, and 
a superffuid bulk viscosity (^3. (The bulk viscosities conventionally called Ci aiid (2 are 
eliminated by the assumption of conformal invariance.) With an additional assumption on 
the form of the entropy current, one can argue that 77^ = 0. Thus, we are left with the three 
transport coefficients 77, k, and (3. 

A holographic superffuid is a superffuid with a dual higher dimensional gravitational 
description. In gravity, and in the simplest possible setup, such conffgurations involve a 
charged scalar ffeld and a f/(l) gauge ffeld, both coupled to the metric [351 Elo Various 
aspects of these relativistic superffuid conffgurations were studied in [6l [71 |39l HQl HH HJl B3l 

mmsisnisz]. 

In this work we construct a solution to the equations of motion that follow from the 
classical gravity system dual to a space-time dependent superffuid and use it to compute 
the transport coefficients described above. The challenge in extending the work of [1] to 
the holographic superffuids of [5] is that [1] relies on analytic techniques whereas the holo- 
graphic superffuid is usually constructed numerically. To simplify our analysis, we work in a 
configuration in which the charge of the scalar, q, is large. We go beyond the probe limit of 
[5] , allowing the metric to backreact on the matter fields. Within our approximation we can 
establish analytically that the shear viscosity of the superffuid component 77^ vanishes, the 
diffusion coefficient k is continuous across the phase transition and that close to the phase 
transition the bulk viscosity (3 scales like the inverse order parameter (and hence diverges). 
For the particular analytic solution of [48] we are able to obtain explicit expressions for the 
third bulk viscosity and the diffusion coefficient near the phase transition which occurs at 



^More elaborate configurations dual to p-wave and d-wave superfluids can be found in [5S1 1571155] . 
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/i/T = 27r/g: 

/ 13 /g/i \-i 21821 - 37152 ln2 
~ T3 \2M^ \T~ ^) ^ 493927r3 

where is associated with Newton's constant in AdS^. It is usually related to the rank of 
the gauge group Nc through a relation of the form Kg oc A''"^. The peculiar dependence of 
Cs on kI will be discussed in section [61 

The remainder of this work is organized as follows. In section |2] we derive the constitutive 
relations for relativistic superfluids, generalizing the work of [2S1 [2H1 EDI 132] • Section [3] 
provides an overview of the computations which appear in sections H] and [51 In section [H 
we discuss the large charge holographic superfluid. The hydrodynamic fields discussed in 
section [H are not spacetime dependent and, with some abuse of language, we refer to these 
configurations as static. Most of the material in the latter section is contained in [SI [IS] 
but we alert the interested reader to several new features of the static superfluid which have 
not been discussed elswhere. These include an extension of the solution of [H] to include 
backreaction of the metric and an analytic derivation of the Josephson condition which was 
observed numerically in [IS]. In section [3 we use a gradient expansion to extend the static 
holographic superfluid solution to a dynamical one and use it to show that diverges near the 
phase transition in section [5?T| that k, is continuous in section [5^ and that r]/s = l/Air and 
rjg = in section [531 Equations and together with corrections to the order parameter 
are derived in section 15. 4[ We end with a discussion of our results in section [HI Appendix 
[XI collects some useful formulae used in the latter sections of the paper while appendix [Bl 
provides an alternate derivation of the viscous transport coefficients by computing two-point 
functions from gravity and employing Kubo formulae. Some of the results of section [H are 
gathered in appendix O 

While this work was nearing completion we learned of |19| which has some overlap with 
the content of this paper. 

2 Relativistic superfluids 

We study a conformally invariant, relativistic superfluid in the limit where the relative su- 
perfluid velocity is small and there are no parity breaking terms in the hydrodynamic expan- 
sion. While our main result for the constitutive relations of the energy momentum tensor 



)+0(f-2.) (3) 
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and current can be found in the literature, we feel that the our construction provides a more 
transparent view of the physical assumptions used. In section 12.11 we discuss the thermody- 
namic variables which describe the superfluid. We then use these variables in section 12.21 to 
construct the energy momentum tensor and charged current of a superfluid in the inviscid 
limit. In section [273] we extend our analysis to include viscous corrections. In section [2^ 
we describe how Kubo relations for viscous transport coefficients can be extracted from lin- 
earized hydrodynamics, and in 12.51 we investigate the speed and attenuation of first, second, 
and fourth sound in a relativistic, conformally invariant superfluid. 



2.1 Superfluid thermodynamics 

As discussed in section [1] a relativistic superfluid can be thought of as a fluid charged under a 
spontaneously broken f/(l) symmetry. In the grand canonical ensemble the thermodynamic 
variables characterizing normal, charged fluids are the temperature and chemical potential, 
T and fi, and their conjugate variables entropy and total charge density, s and pt- Once the 
U{1) symmetry is spontaneously broken then the goldstone boson provides for another 
degree of freedom. In equilibrium, the only extra gauge-invariant scalar degree of freedom 
we can construct from (p is d^j_(f)d^(f) |6]|j Thus, the thermodynamic variables used to describe 
the superfluid can be taken to be the magnitude of the gradient of the goldstone boson, the 
chemical potential, and the temperature: d^(f)d^(f), fi and T. As it turns out, the chemical 
potential and d^(f)d^(f) are not independent. To see this we can turn on an external gauge 
fleld A™* whose time component couples to the charged current the same way a chemical 
potential does, 

C^C + A^^'J^". (4) 

Replacing d^cf) with the gauge invariant combination D^cf) = d^cj) — A"^^ we conclude that 
-Do0 = —pu When Aj^"' = 0, the independent degrees of freedom are: T, p and dicpd^cj) = 
d^j,(l)d^(l) + p? where i = 1, . . . , 3 denote spatial directions. 

We deflne the energy density of the system e(s, pt, di(j)d^(j)) as the time-time component 
of the stress tensor in the rest frame of the normal component. The pressure of the system 
is usually deflned as the force per unit area in an equilibrium conflguration. Since the 
superfluid velocity need not vanish in the rest frame of the normal component an equilibrium 



■^Hcrc, Greek indices run from to 3 and are raised and lowered with the Minkowski metric rj^i, = 
(- + ++)• 

^This relation is also called the Josephson condition. We will see it reappear when we discuss the 
hydrodynamics of superfluids. 
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configuration is not well defined. Thus, we define the pressure P through the relation 



e=-P + sT + fip^ (5) 

and therefore we have 

dP = sdT + ptdp + —d {d^(j)d^(j) - . (6) 

In (ini), X is defined to be the variable conjugate to {d^tpY — p?'. We will see shortly that 
it is associated with the superfiuid density. Our conventions here are slightly different from 
the ones used in [281 El [31]. 



2.2 Inviscid superfiuid hydrodynamics 

Consider ideal superfiuid hydrodynamics where the superfiuid is described by a temperature 
T, a chemical potential a velocity field normalized such that u^u^^ = — 1 and the vector 
quantity d^cf) all considered to be slowly varying functions of the space-time coordinates. In 
what follows we will choose to describe the velocity field of the uncondensed phase and 
d^(f) to describe the non- normalized velocity of the condensate. We define the energy density 
e and charge density pt as the time-time component of the stress tensor T^y and the time 
component of the charged current when in the rest frame of the normal component: 

u^'u^Tf^y = e (7a) 
u^J, = -A . (7b) 

It is standard to decompose the energy momentum tensor and current such that 

(8) 

where j^u^ = j'^^u^ = and ir^^u'^ = 0. Subscripted parentheses denote a symmetrized 
quantity 

Amv) = 2(V + ^-m)- (9) 

The vectors and j'^ describe the energy flux and charge flux in the frame respectively. 
At the inviscid level, the only vector orthogonal to the velocity field is P^^dycf) where P^j, is 
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a projection matrix orthogonal to the velocity field 

Pf,u = rj^u + u^Uy . (10) 

Since we will be using the combination '^d^cj) = d^^cj) + u^u^di,(j) often, we find it useful to 
define two quantities 

Ps = -u^ducj) , yUs?^M = Pl^ducf) . (11) 

Note that 

id,<P){dy) = f,^ {n' - 1) . (12) 

Since (9^<^ is proportional to the superfiuid velocity, we interpret n^^ as the relative superfiuid 
velocity. From the discussion in section [211] we expect that /is = /i. We will soon verify this 
expectation. 

Following ([S]) we can write the charged current as 

(13) 

= (pt-ri/i,)n^' + ri9> 

where ri = Ti{fi,T,dfj,(j)P^'^du(j)). If we assume that carries the charge of the normal 
component then we should make the identification 

Ps = n/is (14) 

where pn = Pt~ Ps is the charge carried by the normal phase. Because of length contraction, 
in a relativistic setting there is no invariant notion of charge density. When we specify 
the charge density, we must also specify the frame. The charge density ps defined above is 
the superfiuid density in the rest frame of the normal component. In the frame where the 
condensate is at rest, the charge density is smaller: p^ = psPs/|<90|. 
The most general expression we can write for T^^ is 

T^u = eu^Uu + P^u {P + T2pln^) + 2r3Psn(^ti^) + Ufijn^n^ (15) 

where P was defined in (jS]), and Ta = Ta{fi,T,dfj^(f)P^'^du(f>), a = 2,3,4. The coefficient of 
P^^ is fixed by the requirement that when the relative superfiuid velocity vanishes = 0, 
expression ( !T5|) reduces to the familiar energy momentum tensor of a normal fiuid. 

To determine X (defined in ([6])), the Ta, and ps, consider the entropy current which 
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is assumed to be carried only by the normal component 



J,^ = sn^ (16) 

and is conserved in the absence of viscous terms. Entropy conservation together with ([5]) 
and (Ej) allow one to reduce 

u^d^T'^" + fid^J'' = (17) 
into a constraint on X, the r^, and fig- Explicitly, one finds 

+ u^u^d^insnf") (rg - fx^Ti) 

+ flsTl" ( -d„T3 + T4,du^s + fJ'du — 



which implies 



+ d,{^,,nn ( -r3 + ^ ) - T2^iy^^u 



^2 = U T3 = r4 = — A = ^ . (19j 

/is /ig /Xg 



and /Xs oc /X. In what follows we make the choice 

fis = fi. (20) 

This freedom in choosing the ratio between /is and /i is due to a freedom of the overall 
normalization of the goldstone boson. The relation (l20l) is known as the Josephson condition. 

To summarize, entropy conservation together with the requirement that the superfiuid 
component carries no entropy implies that at the inviscid level 

T^, = (e + P)u^u, + Pt]^, + 2psn(^P^)"9,0 + ^P/P,"9^09„0 (21a) 

/i 

= -/i . (21c) 
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The thermodynamic relations are given by 



e + P = sT + ^pt 

dP = sdT + pndp — ^d {d„(f)d'^ (j)) ^22) 



2/i 

Pi 



2^ 



= sdT + ptdji — —d{dy(f)d^(t) + fi 

^ Li 

2.3 Viscous superfluid hydrodynamics 

With the definitions ([7]) and the inviscid solution ( 12T|) we can parametrize the viscous cor- 
rections to the energy momentum tensor, current and Josephson condition by V^, r^^, 
and V, 

T^u = (e + PKu, + Pr/^. + 2p,n(^P^)"a„0 + ^P/P^"9;305a0 + 2\/(>,) + ^ (23a) 
J^ = p,u^ + ^P^^d^<P + T^ (23b) 
u^'^^(i) = -{p + v) (23c) 

such that V^j^u^ = T ^u^ = and Tfj^uU'^ = 0. The energy fiux associated with T^j, is 

^ ^ ^ (24) 
= + Ps (<9/,0 - (/i + v)u^) 



and the charge fiux is 



= - id,<l> - (/i + v)u,) + T, . (25) 
p 



There is some freedom in choosing and V^. We can specify the energy fiux seen in the 
rest frame of the normal component (the Landau frame) or the charge density in the rest 
frame of the normal component (the Eckart frame). This amounts to choosing the direction 
of the velocity field relative to or T^. In what follows we define such that 

V; = 0. (26) 

We will refer to this frame as the Landau frame. As in section 12.21 we will use 

p, = -u-'d,^ (27) 
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and 

li,n^ = Pp,cp. (28) 

The viscous corrections to the energy momentum tensor r^y, charge current T^, and 
Josephson condition v vanish when gradients of the hydrodynamic variables can be neglected. 
Thus, it is natural to expand r^^,, and v in gradients of u^, d^(f), fi and T. In what follows 
we will consider these viscous corrections to linear order in gradients. There are many 
distinct expressions which one can construct from the hydrodynamic variables. In principle 
all of these can contribute to r^jy, and v. By distinct terms we mean ones that differ 
after the conservation equations d^T^^ = and d^J^ = have been implemented. In order 
to simplify our analysis of the viscous corrections we will only consider expressions which 
will play a role in a conformally invariant theory and only expressions for which n'^ is small, 

<^ 1. The latter restriction will significantly reduce the number of possible terms that 
can be written down. It is a physically sensible restriction since superfluidity breaks down 
at large superfiuid velocities. 

Within this approximation, the only possible Weyl covariant traceless symmetric tensors 
orthogonal to the velocity field that could contribute to t^u are 

a;, = 2d^,u,^ a% = 29(^n,) . (29) 

Here triangular brackets indicate a symmetric traceless projection onto the space orthogonal 
to the velocity field, 

Kv) = (^A. + ^.a) - , (30) 

and P^^ denotes a projection onto the space orthogonal to the normal velocity as in fllOp . 
Assuming only parity even terms, the only possible Weyl covariant contribution to is 

v; = p,ud''^. (31) 

The only Weyl covariant scalar contribution to v is 

= d,ip,n^) . (32) 

(While not directly obvious, both uf^d^fi/T and d^{psu'^) are equivalent to once we require 
that ^ 1 and that the energy momentum tensor and charge current are conserved.) 



10 



Within our approximation, the most general viscous corrections we can write are: 



T, = -^V; = -KP.^d'^i^ (33) 

One can attempt an alternate analysis of the possible viscous corrections by considering 
the combination 

l^d^J^ + u^d^T'^" = (34) 

and interpreting it as a statement regarding the entropy current similar to what was done 
in ([TB]) - ffTUl) . In the limit ^ 1, it is possible to reorganize the terms of fl25]) in a way that 
makes this entropy calculation simpler. We can redefine Tf^i, and so that in the Landau 
frame equation ( 123|) takes the form 

T^u = (e + P)ui,Uu + Pr]^,u + 2ps/isTi(^Mi.) + Ps/isn^n,, + T^^ (35a) 
Jt^ = ptUfM + Psn^, + (35b) 
-Us = u^d^cj) = -{n + v). (35c) 

In the limit n'' < 1, T^^ ^ and t'"' ^ T^". 

The authors of [211 [22] , found that the form together with the relations 

dP = sdT + ptrfp - ■^d{nln'^) (36) 

and 

duiPsn/,) - d^insTiu) = -d^insUf,) + d^{nsU^) (37) 
(which follows from f l2S]) ) implies 



/i 



(38) 



Following Landau [SO] , if we interpret the term in the parentheses on the first line of f l38p as 
the entropy current and relax the assumption of conformal invariance then we find that 
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the second law, d^J'^ > 0, implies 

= 

V = -Cs^a iPsnn - C2d.u'' (39) 

where 77 > 0, k > and either ^1 > 0, C2 < CiCs or (i = (2 = ^ and Cs > 0. Imposing 
conformal invariance imphes that we must choose the latter which coincides with (l33l) with 
r/s = 0. The observation that (^3 does not necessarily vanish in a conformal theory was made 
in the non-relativistic case in [51]. Before ending this section we note that one should be 
mindful of the ad hoc interpretation of the expression in the parentheses on the first line of 
fl55]) as the entropy current. In [13, [121 dE] it was shown how such an interpretation fails 
in parity violating theories. Also, (136|) is somewhat non standard since it implies that the 
dependence of the thermodynamic quantities like pressure and energy on fi, T and ((90)^ 
gets corrected by viscous effects. 

To summarize, we expect that for <^ 1 

Tf^u = (e + P)UfMUu + P'n^^l^ + ^p^fisU^f^u^) + PsPsn^n^ - w'/.u - Vs(^% (40a) 
= Ptu^, + psUf, - kV^ (40b) 
-/is = u^df,^ = -p + Css' (40c) 

where /ign^ = P'^^dyCp. In the constitutive relations (l40l) . we have kept the full dependence 
on at zeroth order in the gradient expansion. At first order in the gradient expansion, we 
have discarded any terms which contain an n^^ that is not acted on by a derivative. 

2.4 Kubo Formulas 

In this section we deduce some Kubo relations for the retarded Green's functions, 

G^j^'^''{u,k) = i j rf^xe-''="^(t)([T^"(x),r^'^(0)]) 
G''/{u,k) = % j d^xe-'^^'Oit) {[J^'{x),r{Q)]) (41) 
G^^(a;, fc) = I ! d'xe-^^^d{t) ([0(x), 0(0)]) , 
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from linearized hydrodynamics. A similar analysis with somewhat different notation and 
different methods can be found in |3Uj . 

Let us look at small fluctuations about an equilibrium state at fixed temperature, chemical 
potential and zero normal and superfluid velocities, i.e. we write 

T = To + T' /i = /io + /i' = (1, u') = (0, (5>)//i - u') (42) 

where T', fi', Ui and (9j</) are small. At linear order in the fluctuations, the stress tensor (140 ap . 
current fl40bl) . and phase (I40cl) have the form 

Too = e (43a) 
Toi = -(/xpn + sT)ui - psdicj) (43b) 

Tij = (P + '^dku'')6ij - 27]diiUj) (43c) 
Jo = -Pt (43d) 
J^ = u^Pn + - ^ - |9,t) (43e) 

9o0=-/i + C3Ps(^-9,w^) . (43f) 

From f l43cp . we can deduce the standard Kubo relation for the viscosity 77: a velocity 
Ux{y) can be thought of as a small Galilean boost x ^ x — Ux{y)t, which in turn leads to a 
metric fluctuation 6gxy such that dodgxy = dyUx- Thus 

lim - Im G^^'^3'(w, 0) = r] . (44) 

A Kubo formula for k can be obtained by the identification difi = —diAf"^ in f l43ep where 
A""* is an external gauge field. The long wave-length limit of the current-current correlation 
function is then: 

lim lim y Im G'^'Uu, k) = ^ . (45) 

a;-s-0 fc->-0 k T 

One more Kubo relation, also discussed in |10], can be extracted from fl43ep by replacing 9j0 
with the gauge invariant combination di(j) — Af"": 

limReG"'"(a;,0) = . (46) 
Finally, we derive a Kubo formula for ^3. Using current conservation, we write (I43fp in the 
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more suggestive form 



do<P = -/io - - CsdoJ' + ... (47) 

opt 

where by ... we mean second order terms in a gradient expansion. If we turn on an external 
gauge field Aq"' we expect that the left hand side of ( H71) will be modified to do(f) ~ ^o"'- 
can think of the variable as being canonically conjugate to — 9o0. This conjugacy allows 
us to deduce the relations 

\im^ReG'^'^(uj,0) = (48) 
tj-^o tu-^ op 

lim -ImG"^'^(w,0) = Cs • (49) 
2.5 Sound attenuation 

Superfluids admit several types of sound modes [501 ES]. In addition to carrying sound by 
pressure waves, as normal fluids do, superfluids allow for a second sound mode through 
entropy waves. A third sound mode can be generated by surface waves on a thin film of 
superfiuid and a fourth sound mode exists when the normal component of the velocity field is 
prevented from moving. First, second and fourth sound were studied in a relativistic setting 
in the context of AdS/CFT in [HI [391 SO]- In what follows we will study the attenuation of 
first, second and fourth sound due to n, (3 and 77. 

At linearized order in the fluctuations described in f H2|) . the conservation of the stress 
tensor, d^T^^ = and d^T^^ = 0, and of the current, d^J^ = 0, imply that 

{^J^ + ^1^) 9op' + (p^ + doT' + (/ip„ + sT)d,Ui + p,d^<P = (50a) 

Ptd'p' + sd'T' + (/ipn + sT)dQu' + ps9o(9> - ^d'djU^ - 7]dW = (50b) 

3 



^dop' + ^doT' + p^d^u^ + ^9^0 - ^ [d^p' - ^d^r) = . (50c) 

Another important relation in constructing the dispersion relation for the sound modes is 
the derivative di of fl43fl) : 

dodi(P + d,p' - CsPs (^^^ - d.djU^^ = . (50d) 
Assuming the space-time dependence of the fluctuations X'^ = {p',T',u^,d^(f)) takes the 
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form e *'^*+*^^^ we can write flSOap . ( ]50bp . flSOcp . and ( ]50dp schematically as a linear system: 

MX = , (51) 

where M is a 4x4 matrix that depends on u and k. 

As long as the momentum k is small enough, the dispersion relations for first and second 
sound can be determined from the four roots 

U = ±Cak - ik'^Ta + ..., (52) 

a = 1,2 of the determinant of M. To find the roots, we introduce the entropy per particle 
a = s/p and take advantage of conformal invariance. Using the scaling form for the pressure, 
P = T*f{ii/T), we express all the susceptibilities dp^jBT, dpx,/dji, ds/dT, and ds/dfi in 
terms of da /ST. For instance 

3s - PtT— . (53) 



dT dfi fxpt + sT V dT 
At leading order in u and k, the determinant of M has the form 

det(M) = {e - 3co^) (^Ppy - (^p„ + sT)^Lo^^ + ... , (54) 

from which we can read off the sound speeds 

c^ = - c^= (55) 
' 3 ' ' (/ip, + sT) § ■ ^ ^ 

See [IQ] for details. 

The viscous corrections discussed in section 12.31 allow us to determine the sound attenu- 
ation Ta, a = 1,2. By considering the first subleading corrections to det(M) we find that 

Ti = (56) 

for first sound, which is the standard expression for sound attenuation in a conformally 
invariant theory. For second sound we find 

T. _ 2pps ^ fi{m + sT) p,{fipt + sT) 



3{fipt + sT){fip^ + sT) ' 2{da/dT)plT^ 2/i(/ip, + sT) 
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The dispersion relation ( 152|) is valid only as long as the momenta are small Tak <^ c^. 
When ps = 0, the speed of second sound vanishes, and the roots (13^ reduce to 

_u -/ar -7 2 P(PPt + , . 

c = ±c, - zA: r, a; = -^k Jg^rjgfj^s- ' (58) 

provided FiA; ^ Ci. (In this limit, the fourth root sits at u = 0.) 

Fourth sound is associated with wave propagation of the superfluid component when the 
normal component is forced to stay motionless. Experimentally, such a configuration is real- 
ized by channeling the superfluid through a tube packed with fine powder which immobilizes 
the normal component. In this setup the energy momentum tensor of the superfluid is not 
conserved since momentum is transferred to the medium which holds the normal component 
in place. Thus, we consider linear perturbations of the form X = {fi',d^(j)), keeping T and 

= (1, 0) fixed. In this limit, we only need to consider flSOcp and (I50dp which reduce to 

V + ^9^0 - ^d^f^' = 

We assume fi' and di(f) both have a e~*'^*+''^^ dependence. The determinant of the system of 
equations fl59|l is the polynomial 

det(M) =uj^- clk'^ + 2ik'^T^uj + ... (60) 

where 

'' = M^a7^' ^'^me^f^i^' ^^^^ 

and . . . denote subleading terms in the momenta. 

Expanding out the roots for small momenta, T^k <^ C4, one finds 

oj = ±Cik - iTik^ + 0{k^) . (62) 

On the other hand, when ps = we have 

u) = —ik'^—— — TTT^ . (63) 

As expected, these dispersion relations agree with those associated with second sound in the 
limit where sT ^ ppf 
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3 The method of computation 



Holographic superconductors were initially constructed in [5] following the pioneering work 
of [35]. The simplest construction of a superfiuid involves a charged scalar, a U{1) field and 
a graviton, with an action 

'^"2^/^(^^^~ ^^'""^mn - \dm^P - iqAmM'' " ^(IV^H) d^'x + Sb . (64) 

Roman indices m, n = 0, . . . 3;^5 are raised and lowered with the bulk metric. The index 5 
refers to the radial coordinatejj In (164|) gmn is the metric, R is the Ricci scalar, F = dA is 
the field strength of the U{1) gauge field, and ip is a. complex scalar with charge q. The AdS 
radius is denoted L and we will set it to 1 in what follows. 5*;, denotes boundary terms which 
do not affect the equations of motion but will affect gauge theory correlators. At low enough 
temperatures and large enough charge, solutions to the equations of motion following from 
f lM|) admit a condensed configuration where is non vanishing. This is the holographic dual 
of the superfiuid phase. Some introductory material to holographic superconductors can be 
found in [5311511 Eg Eg. 

The goal of the remainder of this paper is twofold: To construct a holographic dual for 
a space-time dependent superfiuid, similar in spirit to the construction of [1], and also to 
compute the transport coefficients r], r]s, k and (3 associated with it. In what follows we will 
summarize our method of computation. 

In the simplest setup, one considers a superfiuid in which the normal component and the 
superfiuid component are motionless. The solution to the equations of motion in this case 
involve the metric (7^1/, Aq, and ip (and perhaps A5, depending on our gauge choice). We find 
it useful to work with gauge invariant variables Gm = Am — dm'P where (f is the phase of ip, 

ij = -^pe^'^^ . (65) 

The AdS / CFT correspondence relates (7^1/ to the metric and energy momentum tensor of the 
boundary theory and Gq to the chemical potential and charge density. A particularly useful 
relation that we derive is that, in our conventions, the boundary value of Gq gives us the 
Josephson condition f l40cp . All these relations will be made more precise in sections HI and El 
and in appendix \M 

^We also use the boundary indices = 0, . . . , 3 and i,j = 1, 2, 3 from before. The indices /i, and i, j 
are raised and lowered with the Minkowski metric rj^^i, and with dij respectively. 
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To obtain dynamics, we need that the boundary theory include non-zero superfiuid and 
normal component velocities. At first we will be interested in an inviscid superfiuid, i.e. 
one where gradient corrections can be neglected or are non-existent. Viscous effects are 
absent if both the superfiuid and normal component velocities are constant in space and 
time. By turning on a constant boundary value of Gj, we introduce a non vanishing but 
constant superfiuid velocity [B]. In the condensed phase, a non zero superfiuid velocity 
implies a non vanishing value for the spatial component of the charged current. Following 
our analysis of superfluids in section [2] we will only be interested in configurations with a 
small superfiuid velocity. Thus we consider only linear perturbations of Gi around the static 
solution described earlier. 

To obtain a solution where both the superfiuid component and the normal component 
are in motion, we Lorentz boost the boundary theory energy momentum tensor and current. 
After such a boost the normal component and the superfiuid component are in motion. Such 
a boost can be achieved in the gravitational dual by a coordinate transformation. 

So far we have a gravitational description of a superfiuid with arbitrary constant super- 
fluid velocity and arbitrary but constant velocity for the normal component. These configu- 
rations have been studied in the literature inP[71[39lllQliIlll2lll3lllllll5lll6lll7j. To go 
beyond these stationary solutions, we need to allow for the velocity fields and other thermo- 
dynamic quantities to be spacetime dependent. When we discussed viscous corrections in the 
boundary theory in section 12.31 we constructed these corrections to first order in gradients 
of the thermodynamic variables. It is natural to carry out the same kind of analysis in the 
gravitational dual. The stationary configuration is described by the metric Qmn, the gauge 
invariant combination Gm and the modulus of the scalar field p. These fields encode in them 
the boundary theory velocity fields and u^, the chemical potential fi and the boundary 
temperature T which are all constants independent of the boundary coordinates x^. If we 
promote all these parameters to be spacetime dependent = n^{x^), = u^{x^), etc., 
then the modified fields gmn-, and p will no longer solve the Einstein equations and mat- 
ter equations. To correct this, we add corrections to the metric 5 gmn, gauge field 5Gm, and 
scalar 5p such that the combinations gmn + ^gmn, Gm + ^Gm and p + Sp solve the equations 
of motion. It is difficult to find the delta'd quantities in full generality, but if we focus only 
on first order gradient corrections to the superfiuid, then it is sufficient to find 6 gmn, SGm 
and 6p which will solve the equations of motion to linear order in the gradients of n^, m^, T 
and p. 

Carrying out this computation we achieve our first goal, to construct a holographic 
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dynamical superfluid. To compute the various transport coefficients, we use the AdS/CFT 
dictionary, described in appendix |X1 The energy momentum tensor of the superfluid is dual 
to Qmn + Sgmn- From it, we can read off the shear viscosities r] and 77^. The charged current 
is dual to Gm and from it we can compute the diffusion coefficient k. Finally, we obtain 
by evaluating the boundary value of u^6G^. 

The computation described above is feasible in principle but in practice technically diffi- 
cult. While we are able to make some general statements about the holographic superfluid 
(e.g. r]s = and diverges close to the phase transition), it is difficult to obtain explicit 
expressions for ^3 and the other transport coefficients. The main obstacle is the absence of 
analytic control over the static solution. To overcome this problem we use the recently dis- 
cussed analytic holographic superfluid of |18]. In [38] an analytic solution to the equations 
of motion were obtained in the limit where the metric does not backreact on the matter 
fields and for temperatures close to the phase transition. In other words, the solution of 



involves a double expansion. One parameter — the charge of the scalar field — allows for 
the matter fields to be weak so that they do not interact strongly with gravity^ The other 
parameter is the distance from the phase transition, T/fi — (T//i)o, with (T//i)o the critical 
value of T//i at which the phase transition occursB Thus, in our formulation there is a triple 
expansion involving the charge of the scalar, q, the relative temperature, T/fi — (T//i)o, and 
gradients, S. 

First, neglecting gradient corrections, in section H] we describe the static solution. To 
leading order in the charge of the scalar, the metric completely decouples from the matter 
fields and the solution to the Einstein equations is the Schwarzschild black hole which we 
describe in 14. 1[ At the next order in q, we solve for the matter fields. This is described in 
section 14.21 We then go beyond the probe limit in section 14.31 where we consider the leading 



order backreaction of the metric. Up to section 14. 4[ our analysis is general. We express the 
particular analytic solution of [IH] as an expansion in T//i — (r//i)o in section 1^31 Section 
[5] extends the solution of section H] to include gradient corrections a la [Ij. The extension is 
performed order by order in q in sections |5. II to [531 The principal results are that ^3 diverges 
at (T/yu)o, K is continuous across the phase transition, rj/s = l/47r, and rj^ = 0. For the 
analytic solution of [IS] , the explicit values of and k described in ([2]) and are obtained 
in a perturbative expansion near T/fi — {T/ fi)o in 15.41 



^This same "probe limit" was used in one of the first holographic superconductor papers [5]. 
''In [IS], the magnitude of the order parameter replaced T//i — (T//z)o as the small parameter of the 
expansion. 
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4 Static holographic superfluids 



Before describing the solution to the equations of motion that follow from (!64|) . we need to 
explain the precise form of the boundary action Sb and also the relationship between the 
bulk gravity fields gfj^u, and ip and the boundary field theory quantities T^" , J^, and (p. 

The boundary action Si, can be obtained by requiring that the variational principle is 
well defined. Assuming that the metric Qmndx'^dx"' becomes asymptotically anti-de Sitter 
at large r, 

lim ds'^ = -r'^dt^ + S^^dx'f + 2drdt , (66) 

i 

we take Sb to live on a constant and large r surface. On this surface, we can define an 
induced metric 7mn and a unit vector (the lapse function, not to be confused with the 
relative superfluid velocity n^) that points toward larger r. The terms in Sb relevant for this 
paper are 



1 



2k? , 

5 boundary 



7 {2K - 6) d'^x 

(67) 



+ i /" ( ^aI^I' + \rhl {Vn^dm^ + c.c) ) d^x . 



5 boundary \ 

Here K^J,y is the extrinsic curvature, Ky,y = — F^^/ \/g^ with the Christoffel symbol, and 
K = K^,-i^^\ 

The expressions for m/^ and ttia depend on the scalar potential 

Vm')=m'\^\' + Om'), (68) 

in ( l6l|) . Through the AdS/CFT dictionary, the mass of the scalar is related to the con- 
formal dimension A of the dual operator through 

= A(4 - A) . (69) 

By requiring a well posed variational problem and also that the boundary action is finite, 
one finds that and rfiA in (l^7|) must satisfy 

and »iH . (TO) 

2 A < 2 
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The prescription for obtaining the energy momentum tensor in the dual (boundary) gauge 
theory, T^^, was derived in \57\ following the earlier work of [21 E]- For configurations in 
which the fields depend only on the radial coordinate r, 

i^l{Tf,u) = - lim | K^^ - K-f^,^ + 87^^ - ^t^j^'^aIV'P 

- ^If^urhA {iJ*n^'df,tp - 2tp*n(^f,d^)ip + c.c.) ^ . (71) 

To understand the relationship between and ip, note that ip has the large r expansion 
^ = V'.r-^^-^) (1 + . . .) + Ar-"^ (! + •••), (72) 

with . . . denoting subleading powers of The coefficient ips is associated with a source 
term for the operator and the coefficient ipr is associated with (O^). If we set ^/^^ = 
implying that we are not deforming the boundary theory, and we consider configurations 
which depend only on the radial coordinate r, then the expectation value of the operator 
dual to ip is given by 

= v^Ca^. , (73) 

where Ca denotes a dimensionless constant whose value depends on our conventions for 
normalizing the source term ipg [59]. Its explicit value will not play an important role in this 
work. We split ip into its modulus p and phase (f as in ( 165|) . Following (172|) and (!73|) . we 
identify the goldstone boson with the phase of (in the absence of sources) through 

(0)=lim<^, (74) 

r— >oo 

and the modulus of (in the absence of sources) with the near boundary asymptotics of p, 

kI\{0^)\ = Ca lim r'^p. (75) 

r— >-oo 

The boundary theory current J^^ is related to the bulk gauge field through 

2/T\ T 2 A ( possible \ 

H^{Ja) = lim r y4„ + divergent 1 (76) 
1 S-OO V terms / 



^For certain values of A the series expansion associated with tps and the series expansion associated with 
ipr overlap. In this case one obtains logarithmic terms in the series expansion. We will see an explicit 



example of this sort of behavior in section 14.51 
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which is vahd for gauge fields which depend only on the radial coordinate r and the gauge 
choice = 0. If the configuration admits an event horizon at some rb = 1 then the chemical 
potential of the boundary theory can be obtained from 



The results ( I7T|) . ( I73|) and ( 176!) can be modified to include more intricate spacetime 
dependent configurations by including appropriate boundary counterterms. In general, we 
expect these counterterms to depend on two or more derivatives of the bulk fields in the 
directions transverse to r. For the analysis carried out in this work, all such terms can be 
neglected and we will use f lTT]) . fl73|l and fl76|) from now on. 

Bulk configurations with ip = and Aq ^ correspond to boundary theory configurations 
with a non vanishing charge density. Once is non vanishing and regular, and ^/^^ = then 
the boundary configuration corresponds to a spontaneously broken phase of the theory where 
0, the phase of O^, is the goldstone boson. 

To simplify our computations we formally expand the metric and matter fields in inverse 
powers of the scalar charge and work in the large charge {q oo) limit: 



9.rn = 9^1 + g^lq-' + 0{q-') A^ = Agk~' + 0{q-') = + 0{q-') . (78) 



The ansatz fl78p is useful since the bulk energy momentum tensor of the matter fields is 
quadratic in ip and A^, effectively decoupling the Einstein equations from matter at leading 
order. In sections 14. ![ 14.21 and 14.31 we will discuss the construction of a superfluid order 
by order in the large charge expansion. In section 14.41 we extend the solution to include 
a non-trivial velocity for the normal component. In section 1^31 we show how to construct 
an explicit solution to the field equations to order close to the critical temperature. 
We emphasize that in this section all viscous corrections to the stress tensor and current 
vanish because the solutions presented correspond to a fiuid moving with uniform superfiuid 
velocity, normal velocity, temperature and chemical potential. 




i/b 



(77) 
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4.1 The AdS-Schwarzschild black hole. 



At order g° the solution to the Einstein equations is the AdS-Schwarzschild black hole given 
by the line element 

ds^ = r^-{l-^)dt^ + y^idx^ ] + 2drdt . (79) 



- dt^ + J2{dx'f^ + 2drdt . 



The black hole horizon is located at r = 1/6 and the asymptotically AdS boundary is located 
at r — oo. The Hawking temperature of the black hole, which is also the temperature of 
the boundary theory is given by 

We have chosen to write the metric in ingoing Eddington-Finkelstein coordinates in order to 
avoid a coordinate singularity at the event horizon. 

Using the prescription f lTT]) we find that the configuration f l75]) corresponds to an energy 
momentum tensor 

(T,,,) = diagonal (e |e ie |e) (81) 

where 



with b given in 



4.2 The probe approximation 

At order q~^, the equations of motion reduce to the Maxwell and Klein-Gordon equation 
for ip and Am- We start with an ansatz where A^"^ and A^^"^ are non zero, and the spatial 
components of the gauge field Ai are turned on at the linear level. Consider the gauge 
invariant variables]^ 

Gm = Am- dmf ■ (83) 

^Since we are working to leading order in q in the matter fields, we can use Am and ip instead of 
and t/j^^y Alternately, we can omit the superscript (1) from all expressions and reinsert them only when 
discussing boundary observables. 
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If we keep only the term in the scalar potential (!68|) 



we obtain the equations of motion 



dr \ dr J b^r 
d^p hh'^r^ -1 dp _ b'^m'^r^ _ _ b^r^p Gl 

+ -l)lt~ b^r* ~{b^r^ - 1)2 ' ^ ' 



which are supplemented by the constraint 

In the more conventional Fefferman-Graham coordinate system it is customary to choose a 
gauge where the scalar field is real in which case one also finds that = 0. From (!83|) and 
(!85|) it is clear that in the Eddington-Finkelstein coordinate system one can not set both 
if = and = unless the scalar field vanishes entirely. Of course, once the solution in 
the Fefferman-Graham coordinate system is known, one can obtain the relation between G^ 
and Go appearing in fl85l) by an appropriate coordinate transformation. 

The boundary conditions we impose on our fields are that Gm and p are finite at the 
black hole horizon. Near the asymptotically AdS boundary we require that the source term 
for the scalar field ips, defined in (1721) . vanish. From the constraint equation ( l85l) we see that 
finiteness of G5 at the horizon implies that Go must vanish there and finiteness of G5 at the 
boundary follows from finiteness of Go. Thus, 



p{l/b) = finite Go(l/6) = , (86) 

and 

p > 0{r-^) lim Go = finite . (87) 

s>oo r— >oo 

It is important to point out that we did not require, a priori, that Gq vanish at the 
horizon. Rather, this restriction followed from the constraint equation (185|) . When we 
consider a time- dependant (dynamical) superfiuid in section [5] then the constraint equation 
equivalent to (1851) will imply that Gq is non vanishing at the horizon. When working in 
the Fefferman-Graham coordinate system it is standard practice to require that Go (or Aq) 
vanish at the horizon on the grounds that Godt have finite norm there. When working in 

^"^In principle, one can carry out an appropriate rescaling of the couplings appearing in the scalar potential 
so that higher orders of tp will also contribute to the equations of motion at order q~^. This has been done 
in, for example, [39j. 
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ingoing Eddington-Finkelstein coordinates, a finite value of Gq at r = 1/6 will not generate 
a divergence on the future horizon and so, if there are no other constraints, we are free to 
keep its horizon value arbitrary. The norm of G^dt does diverge at the bifurcation point of 
the horizon if G'o(l/&) 7^ 0. However, when constructing gradient corrections to the black 
hole geometry, the past horizon (which includes the bifurcation point) becomes singular so 
having Go diverge there is not a cause for worry. The singular nature of the past horizon is 
not surprising since generic solutions of viscous fluid dynamics are not expected to be regular 
in the infinite past. 

Once a solution to (184|) is obtained it is convenient to choose a gauge where = 0. In 
this gauge, the phase of the scalar ip is given by 

dr(p = -G5 df,ip = - lim G^ (88) 

r— s-oo 

where the second term in f l88p comes from the requirement that is not sourced, by which 
we mean limr_!.oo = 0. Our definition of if allows for an arbitrary additive constant which 
has no physical significance. With the gauge choice fl88|) . the bulk to boundary identifications 
(EH) and (I76D take the form 

-9i(0)=limGo (89) 



and 



4{J.) = hm r^G. + ( i:™) . (90) 

r— >00 \ terms / 



Using the definition of the chemical potential f l77p . our solution for (p and the horizon 
boundary condition Go(l/&) = 0, we find that (18 9 p implies that 

dM = -/i (91) 

which is precisely the Josephson condition fl40cp in the boundary theory in a configuration 
where the spatial components of the normal velocity vanish. Put differently, by requiring 
that Go and G5 are finite at the future horizon and that Aq is not sourced, the horizon 
asymptotics of the constraint equation fISS]) enforce the Josephson condition in the boundary 
theory!"] 

By turning on Gj we extend the solution to include a non vanishing superfluid flow 
di{(j)) 7^ in the boundary theory. Working with a small value of Gj corresponds in the 
current setup to setting n'^ ^ 1 in the boundary theory. As discussed in section m this 



'^A numerical derivation of the Josephson condition when q is finite can be found in |45j . 



25 



limit is physically sensible since superfluidity breaks down at large superfluid velocities. The 
linearized equations for Gi are given by 

dr^ rih^r^ - 1) dr r%^ -I ^ ' 

Multiplying equation fl92|) by Go and using f l84ap . the equation for Gi takes the form 

d fr%^-l ( dG, ^ dGo\\ _ G.dGo 

<^o-r--<^*-r- ---T-j—- \^'^) 



dr \ Ar \ dr dr J J dr 

We will find that the form (1931) is more useful than ( l92l) when discussing the backreaction of 
the metric in section 14.31 At the horizon we require that Gi is finite. To understand what 
boundary conditions to impose on Gi at the asymptotically AdS boundary we note that in 
a gauge where = 0, (IHHjl implies that 



di{(P) = - \im Gi . (94) 

r—>oo 

Since (p2l) is linear, it is convenient to define g{r) such that 

G, = -g{r)dM (95) 

where is a constant and g satisfies the same equation as Gi but with boundary condi- 
tions such that g{oo) = 1. In [51 [71 HSl HHl SS] the equations of motion for Gi were studied 
beyond the linearized approximation. 

4.3 The backreaction of the metric 

At order the metric gets corrected due to the matter fields. The most general isotropic 
order gauge-fixed line element is 

g^^dx'^dx'' = -r^f^^ {r)dt^ + 2s^^^ {r)drdt . (96) 

The equations of motion for f^"^^ and s'-^^ are 



dr (6^4 -1)2 \dr^ 
dr^ ' \ dr J ¥ \dr J ^ 
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We require that s^^^ vanish at the asymptotically AdS boundary and that the solution for 
/••^^ does not shift the mass of the Schwarzschild black hole b in (17^ . 



hm rV^2) = + fd':™^o':t) . (99) 

1 ^OO \ terms / 

(2) 

Once we turn on Gi at the linearized level, we must allow for fluctuations of the = 
^'^iT'^ components of the metric [iQ],!^ This coupling represents the physical fact that a 
charged current should carry momentum in a medium with nonzero charge density. The 
linearized Einstein equations for 7^ are 

- ir"'^^ = -r'^^ - ^'^'^-/^'^o . (100) 
dr \ dr J dr dr h'^r^ — 1 

Because the '-yf'^ are coupled linearly to the Gj's, in analogy to ( p5|) it is natural to define 7 
through the relation 

7^^-7 5,(0). (101) 



After a manipulation similar to the one that took us from ( l92l) to (1931) we obtain 

^ = ir'/-^] . (102) 



dr \ dr J dr \ dr ^ 
Using equation (!93l) we find that the solution to fll02p is 

r%^ - 1 ( dg dGo\ K^n /,4 4 ,\ nnq^ 

(The manipulations leading to fllOSp essentially follow from the work of [12].) The integration 
constants in (11031) were fixed as follows. We require that 7 is finite at the horizon and that 
near the asymptotically AdS boundary it satisfies 

T 4 ^ f possible \ /-ir\A\ 

hm r 7 = Ps + I divergent 1 . (104) 

) >-00 2 V terms / 

The boundary condition ( 1104^ follows from the linearized version of the inviscid expression 
i^The 7f' and 7 in this section are unrelated to the boundary metric 7^^/ of section |4l 
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for the energy momentum tensor (12T|) . 

(Toi) = {e + P - fips)uoUi + —do(j)di(j) 

(105) 

and the relation ( !7T|) between the boundary theory energy momentum tensor and the bulk 
metric. The freedom we have in choosing the overall coefficient of the term in a series 
expansion of 7 is due to possible shifts in the spatial components of the normal fluid velocity 
which we must set to zero. 

It is interesting to note that the condition fll04p implies that the g^i components of the 
metric vanish at the horizon: 

7(1/6) = 0. (106) 

If instead of fll04p we had used a different boundary condition for 7 then this would have 
implied, via f ll05p . that we had turned on the spatial component of the normal velocity Ui. 
Once Ui is non zero the spatial component of the entropy current defined in (fT6|) will not 
vanish and we can expect from the analysis of pD] that gQi{l/h) ^ 0. Thus, requiring that 
= implies gQi{l/h) = 0. 



4.4 A boosted static solution 

The solution in sections 14.1114.31 is static in the sense that it is time- independent: the energy 
density and charge density are constant in spacetime and so is the velocity of the superfiuid. 
It is straightforward to extend the static solution described so far to include a spacetime 
independent velocity field for the normal component. Let A^,^ denote a boost by a velocity 
parameter —u'^ where 

""^-y ^ ' 

and = X] &i < 1- Under the coordinate transformation — x'^ = K^^Xp, the field 
transforms into 



G' 




-Gqu^ - gN^ 
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where we defined 

iV^ = A^*9i(0) . (109) 

The boosted metric takes the form 



1 f^^^\ ( s^^^ 

ds' 2 = -r^ I 1 - j^^^ + u^Uucix'^dx" + r'^P^^dx^dx" - 2 j^l + — ) u^dx^dr 

+ 2rV(M^t')^a;^^a;'' • (HO) 

In flllOp P^jy is the projection matrix defined in ffTOj) and s'-^^ and f^'^^ are the order 
corrections to the metric defined in ( l96|) . The scalars G5 and p remain invariant under the 
coordinate transformation. We emphasize that f lllOp and f llOSp solve the equations of motion 
to order since they have been obtained via a coordinate transformation of solutions. 
In terms of boundary theory observables, we have 

= - hm (111) 

so that according to fllOSp 

N^ = - lim GoM^ + 9^,(0) 

= -/i«,+9:,(0) (112) 

where in the second line we used fl9Tl) and in the third line we used the fact that since the 
configuration we are considering is time independent /is = p- In the rest of this work we will 
consider only the boosted frame and will omit primes from boundary quantities which are 
all boosted. In the boosted frame the charged current is given by fl21bp with 



2 V 2r^ , f possible 

KgPt = — hm r Go + I divergent 
r— >00 V terms 

2-1 r 2 , / possibi 

'^S/^ Ps = — hm r g + [ diverge, 
r—>00 V terms 



(113) 



The energy momentum tensor in the boosted frame can be read off the metric flllOp using 
([71]). We find 

{Tf,u) = eUf,u^ + i^(^Piiv + 2psfiu^^n^) (114) 

where e and P^i, were defined in (182|) and (fTOj) . This is the expected form for the inviscid 
energy momentum tensor (1210 when working to linear order in n^. 
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To summarize, our strategy for computing the bulk dual of a static holographic superfluid 
is to solve and (P7|) to obtain the matter fields Gq, p, and the subleading order 
corrections to the metric s*-^-* and /'•^^ With these expressions we can solve the linearized 
equations fl92|) and f llOOp for the spatial component of the matter fields Gi = —gdi{(j)) and 
the space-time component of the metric Goi = —r^jdi{(f)). 



4.5 An explicit solution near Tq 

In [IH] it was shown how to construct an explicit solution to fl84p . f l85l) and fl92l) close to the 
phase transition. In what follows we reproduce the solution of |18] and extend it to include 
the leading order backreaction of the metric. 

Working in the probe limit, we denote the solution to ( 18^ in the normal phase by pG^Q^ 
where according to fl77|) fi is the chemical potential of the boundary theory and 

.(0) ^^ 1 



G'o' =[^-^)- (115) 

Condensation of the scalar field implies that at a given temperature there exists a critical 
value of the chemical potential (which we denote by = /iq) at which the equation of motion 
for p in (I84bl) admits a zero mode. By a zero mode we mean a solution to (I84bl) but with 
Go replaced by poG^q\ This solution, which we denote p'^\ is defined up to an overall 
multiplicative constant. Working perturbatively in the dimensionless quantity (u — po)h it 
is possible to obtain the solution to f lH^ close to the critical chemical potential^ Such an 
expansion takes the form 

oo 

6Go(r) = p^hGf\hr) + Y,Gf^\hr){ph - pobT 

oo __ (116) 



n=l 



The boundary condition Mnir^oo Gq = p implies that lim^^oo Gq"^ = for n > 2 and 
limr_>.oo Gg^^ = 1. We have conveniently defined the argument of p*-"^ and Gq"^ to be the 
dimensionless combination br. 

In practice the expansion (11161) has been found useful only when the zero mode of p can 



^■^In [48] a slightly different expansion was carried out where the small parameter was (O^). The critical 
chemical potential can then be determined as a function of (O^). See appendix |B] for details. 
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be obtained in closed form. Such a closed form solution exists for the special case = — 4jl^ 
where the phase transition occurs at 

/io6 = 2. (117) 

Plugging the expansion flll6p into the equations of motion (l84l) and expanding in powers of 
jj,b — jiob one finds the following equations of motion for Gq"^(x) and p^^\x): 

(x^Gf'(x))' = sT\x), (118) 
+ = Sf, (119) 

where S'q"^ and S'p"''' are functions of the lower order solutions p'-"^) and Gq™'' with m < n. 
The first few solutions for Gq"^ and p*^"^ are given in appendix O 

Once the solution to (184|) has been obtained perturbatively, it is a simple exercise to 
compute the solutions to the linearized vector equations for Gi, ( 192|) . In the notation of (p5|) 
we expand the vector modes near po^ = 2 such that 

oo 

9{r) = Y,g'^^\rb){fib-2r. (120) 

n=0 

Imposing the boundary condition g{oo) = 1 amounts to setting lim^-^^oo fi'^'^'* = 5"°- Inserting 
the expansion (I120p into fl92|) . we find that the equation of motion for g^"^^ takes the form 

(^^^?(^)'(x)) =5f). (121) 

where S^""^ depends on g'''^\ Gq"^^ and p^™-* for m < n. The first few values of g^"'^ can be 
found in appendix O 

Using the explicit solutions for g^^\ G^q''^ and p'-"-' and the bulk to boundary mapping 
described in sectionlHand summarized in appendix|Al we can compute | (O^) | and the bound- 
ary theory current (J^) for a configuration moving with a superfiuid velocity proportional 

^"^The case = —4 corresponds to A = 2 where the series expansion ([7^ reads 

-(/, = i/j.r-^ Inr (1 + . . .) + Ar^'^ (! + ...)■ 



31 



to ^^'{(|)). More explicitly, using ([90]), (JTS]) and ^ with m'^ = (1,0), we find that 

nlh'^p, = 2 + 7(/i6 - 2) - ^(96 ln2 - 71)(/i6 - 2)^ + C - 2)^) 

/ 241 \ 

nlb^fi-^Ps = 6(/i6 - 2) - ( 96 In 2 - — j - 2)2 + C - 2)^) 



(122) 



and 



1^1(0,) I = 4v/3(^6 - 2)V2 + ^^^^^^ ^^^ _ ^ ^ ^^^^ _ ^^^3) 



Similar expressions can be found in |48] . 

The equations of motion for the backreacted metric are similar to those of the matter 
fields. Expanding the corrections to the metric s^'^\ /'•^^ and 7, 



00 

= (2) 



^s(2'2")(r6)(^6-2)" 

00 

/(2) = fi^^) {rh) {fib - 2)" (124) 

n=0 

00 

7 = 6^7(2^)(r6)(/i6-2)" 



n=0 



and inserting them into the equations of motion fl97|) and into fllOSp . we obtain a set of 
equations for s^'^'"'\ /(2'") and 7*^"'^ The solutions for the first few (n) are given in appendix 

o 



5 Dynamical holographic superfluids 

To generate a spacetime dependent holographic superfiuid we can follow the same strategy 
that led us from the inviscid superfiuid in section 12.21 to the viscous superfiuid in section 
12.31 we allow the hydrodynamic variables to depend on the spacetime coordinates and look 
for the appropriate corrections to the bulk fields. This is essentially the strategy used by 
[1] to construct the metric dual to a viscous fluid. In section 15.11 we work out the dual of 
viscous flow to leading order in the large q limit. Since at leading order only the metric 
enters into the equations of motion, this is essentially a rederivation of the results of [1] . The 
reader familiar with the work of [1] may go directly to section 15.21 where we compute the 
effects of viscosity on the bulk matter fields. In section 15.31 we argue that the shear viscosity 
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associated with the superfluid motion vanishes. Finally, in section 15.41 we obtain an explicit 
expression for the third bulk viscosity, diffusion coefficient and viscous corrections to the 
order parameter close to the phase transition using the solution in 14.51 

5.1 Viscous flow from the Schwarzschild black hole 

Consider the solution flllOp to the Einstein equations but with s^^-* = /'•^^ =7 = 0, i.e., 
the boosted AdS-Schwarzschild black hole. The boundary theory stress tensor dual to the 
configuration in fill 01) is given by (11141) with Ps = which describes an inviscid fluid moving 
at uniform velocity u'^. Viscous corrections to the fluid motion vanish since gradients of the 
velocity fleld and energy density are zero. 

To obtain a conflguration where the velocity fleld and energy density are not uniform we 
promote /3j in f ll07p and the inverse temperature b in f lTS]) to be spacetime dependent. In 
doing so f lllOp is no longer a solution to the Einstein equations and the metric has to be 
corrected. We denote the correction of the metric by Sglnn and, following [61], parametrize 
it by: 

^g^Z.dx'^dx'' = -2u^dx^'r + dx" - 26s^^^u^dx>'dr + r^6h^°^ P^^dx^dx" 

- r'^Sf^'^^u^dx^u^dx" + r'^26V^fu^)dx^dx'' + r^STrfJdx^dx'' . (125) 

Here 57r(0) is symmetric and traceless and both SV^^"* and ^vr*^"-' are orthogonal to the velocity 
fleld n^. Circular parentheses were deflned in and denote a symmetric combination. We 
will flx most of our gauge freedom by setting Sh^^^ = as in [62]. At this point it would 
perhaps be useful to recapitulate our notation. In section |4] we used an unbarred superscript 
to denote coefficients in a series expansion in the inverse charge 1/q. In section 1^5] we used a 
barred superscript to denote coefficients in a series expansion in the chemical potential. Now 
we use a 5 to denote gradient corrections to the metric. Our notation differs from the ones 
used in [H |62l [63] where a superscript denoted coefficients in a gradient expansion and from 
the one in [62] where a barred superscript denoted coefficients in a near boundary expansion. 

We will not attempt to solve the Einstein equations entirely. Rather, we neglect all 
terms which involve two derivatives of the velocity fleld and inverse temperature and look 
for a solution linear in these gradients. We will denote the terms that have been neglected 
by Oi^d"^). This means that the metric components in fll25p depend only on one derivative 
of b or u^. The equations of motion for Ss^^\ Sf^'^\ SV^'^^ and 57i^^^ naturally decompose 
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themselves into scalar, tensor and vector equations under the 5*0 (3) C SO {3, 1) under which 
is invariant. They are given by 

= 

dr 

dr ^ 

(126) 





d 


f,d6Vi'^ 
\ dr 




dr 


d 

dr 




dSirj^J 
^ dr 



with 0"^^ defined in (1291) . In addition to the equations of motion (11261) there are four constraint 
equations 

u>'d^b = hd^u^" P^'^d^b = bu^d^u" . (127) 

(Equations (I127P amount to four independent equations since the matrix P^^, is a projection.) 
The constraint equations are equivalent to energy momentum conservation d^T'^'^ = with 
an energy momentum tensor as in (I114p but with ps set to zero. The boundary conditions 
we impose on the metric components are that they are not singular at the horizon which is 
located at r = 1/6 and that they do not deform the boundary theory metric. The requirement 
that u^u'^T^y = e sets the (9(r~^) term of the near boundary series expansion of 5f^^^ to 
zero. Working in the Landau frame ( l26l) sets the 0{r~^) term of 5V^^^ to zero. We refer the 
reader to [U [62] for more details. 
The solution to (I126P is given by 

(128) 



(^-^-^ arctan (rb) - \n{rb) + ^ ln(l + r6) + ^ In (l + r%'^)^ (T^^ . 



Following f[7T|) the boundary stress tensor is given by 

I ,. . 1 



^4{T,u) = ^ i^u.u, + V,.) - (129) 



Using the Bekenstein-Hawking entropy formula we obtain 

2nb^ 



4 



(130) 
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which gives us the celebrated relation ([T]). 

For the reader unfamiliar with the formalism introduced in [1], we point out that a key 
feature of the equations of motion which allows one to simplify them considerably is that 
they are ultra-local: one can solve equations (11261) in a neighborhood of a point Xq and patch 
these solutions together. Technically, this feature of the equations of motion allows one to 
work in the neighborhood of a point Xq where the velocity field and inverse temperature are 
chosen to satisfy 



u 



fl o) + x^^X + (^(-9') 

^ ^ (131) 

b = bo + x>'d^b + 0{d^). 

When working with a computational software program such as Mathematica the decompo- 
sition (I13ip significantly simphfies computations. We refer the reader to |1] for details. 

5.2 Viscous superfluid hydrodynamics in the probe Umit 

To compute the viscous corrections to the current and to the Josephson condition we follow 
the same strategy as the one outlined in section 15. 1[ In addition to promoting b and to 
spacetime dependent quantities, we also allow the chemical potential fi and the field A^^ (of 
which we keep only linear terms) to be spacetime dependent. Eventually we will be interested 
in configurations where is small so we will set A^'* to zero but not its derivatives^ We 
then look for corrections to Gm, p and ip which we denote by 6Gm, Sp and 6Lp. It is convenient 
to decompose 6G^ into terms parallel and orthogonal to the normal velocity u^: 

6G^ = -6Gu^ + 6g^ (132) 

with u'^Sg^^ = 0. We then require that the Maxwell and Klein-Gordon equations of motion 
be satisfied to first order in gradients of the hydrodynamic variables. 

As in the case of the Schwarzschild black hole, we can use the ultralocal nature of the 
equations for 6p, 6G and 5(7^ to solve the equations of motion around a point Xg = where 



^^To be precise, going back to (|112p . A^^ receives viscous corrections because /in^ ^ —fiu^ + d'^{(j)). Since 
is small and we are working to linear order in gradients these corrections do not play a role in our 
computations. 
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the fields X = {Gq, G^, p) can be expanded in the form 



X(r;6(a;^),/i(a;^))=X(r;6(0),/i(0)) + 



9X 
'db 



dx_ 

dp. 



x'd.p + Oid^) . (133) 



x=0 



Because we work in a limit where n'^ is small, the expansion for G* is simpler: 

G,{r;b{x^),p{x^),d,<f){x^)) = -g{r;b{x^),p{x^))d^<f){x^) = -g{r;b{0), p{0))x''d,d^<f) + 0{d^) 

(134) 

The transformation properties of 6G, 6g^, 6G5 and 6p under the 5*0(3) C 50(3, 1) under 
which is invariant imply that the equations of motion for the scalars 6G, 6G^ and p will 
decouple from the equations of motion for the vector 6gf^. 

5.2.1 Scalar sector 

In the scalar sector the equation of motion for 5G and 5p are 

fj4j.5 dr ^ dr ) P ^ P P ~^ r(r^M — 1) " dr ° 

3{r%'^-l)\ r(r464_i) h ^ y u 

(135) 



and 



+ 1 - 2r%^ f d 



158 



,5p 



56V -1 d 



dr"^ ' ' r{b^r^ — 1) dr 
2pG6G - G^6p 



5p 



b^m^r"^ 
b^r^ - 1 



5p 



r%^ - 1 / _ d ^ 

3Vp + 2r—Vp 

dr 



^354 



^1. a L 2 4r^bVO 5r%' + 3 d 

+ -OaU 2m rp — — — — 

3 V ^ b^r^ - 1 r%^ dr 



(136) 



In both (I135P and ( 1136^ we have used 



VX = —u daP + 7;^daU 
op 3 ob 



(137) 
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with X a bulk field. The constraint equation is given by 



r%^6G d.N^dg d 



- 1 Vp^ dr 3(r464 _ 1^ 

In presenting fll35p -( !T38|) we have used (11270 . 

The boundary conditions on 5p and 5G are that they are finite at the horizon and that 
the source term ipg in ( 172]) vanishes. We impose 

hm r^^G = + f i:r"ent) (139) 

r— >-CXD V terms / 

which amounts to the requirement that u^J^ = —pt- To the order we are working in, we can 
also consistently set 

lim {G' + 5G^) = hm {-G^u^ - gN^ + 5G^) = -9^0 . (140) 

The viscous corrections to 9^j<^ are encoded in the boundary value of 5G^. By contracting 
(I140p with we obtain an expression for the viscous corrections to the Josephson condition, 

- u^d^cj) = lim (Go + SG) . (141) 

r— >oo 

Here, as in (llOSp G'^ denotes the value of Gfj_ in the boosted frame. In principle, since in this 
section the velocity of the normal component is space-time dependent and generically does 
not vanish, we could have used primed variables such as 6G'^ instead of the current ones. 
We have not done so in order to avoid cluttering the notation. 

Once we take to be small, the near boundary asymptotic expansion of 6G5 takes the 
form 



mw4 



= :r"' , {ny, + nVt) + 0(r2^"^) . (142) 



1 2 ^,.4 

Finiteness of the action implies that 5G^ should fall off faster than near the boundary. 
One way to see this is to consider a gauge where = 0. In this gauge, the boundary value 
of the variation of the action with respect to would receive a divergent contribution from 
5G5 unless it falls off fast enough at large r. Thus, (I142p together with ( 12T|) imply current 
conservation in the boundary theory. 

The constraint equation f ll38p also supplies us with the correct horizon value of 5G. Near 
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the horizon, the constraint equation fll38p together with regularity of Sip imphes that 



5G{l/h) - ^ 



12 dr 



OyU H — -V 



rb=l 



dr 



(143) 



rb=l 



To obtain f ll43p we expanded Go and p near the horizon located at r6 = 1. As opposed to 
the static configuration where we had GQ{l/h) = 0, here 6G does not vanish at the event 
horizon. Since 6G{l/b) oc p(l/6)~^ one might erroneously conclude that 6G{l/b) diverges in 
the yU — yUo limit where p vanishes. To see that this is not the case consider the solution to 
the equations of motion near the phase transition. As discussed in section 14.51 the solution 
takes the form 

Go = ^i(l- 4t^] +0{pb- flob) (144) 



and 

p = 0{ifib-fioby/') (145) 

for any value of the mass of the scalar field. Inserting (11441) and fll45p into f ll43p we find 
that 

5G{l/b) oc /^o^.^; + 3^;^./^ . (146) 
/i6 - /ioO 

To leading order, current conservation close to the phase transition amounts to 

= (pn^O 

dp , , dp 
= —u dub + —u dup, + pd Uy 
db dpi 

2p p (147) 

b /i 

\d''uy + —u^dyp] p. 
3 po J 

In the third line we used the fact that near the phase transition, p oc b^^ (which follows from 
p oc p and dimensional analysis). In the last line we have used fll27p . Thus, close to the 
phase transition the leading divergent contribution to 6G vanishes and it follows that 

SGil/b) = 0{{pb~pob)') . (148) 

Without an explicit solution to the equations of motion f ll36p and (11350 it is difficult 
to compute the viscous corrections to the Josephson condition and to the expectation value 
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of Ojjj. However, it is possible to argue that "will diverge close to the phase transition like 
the inverse power of the order parameter squared, or {fib — /igfe)"^. Equation f ll35p together 
with fll48p suggest that bG should scale as — fiob)^- Thus, we should find that close to 
the phase transition, 

u^'d^<j) = -/i + O {{fib - fiobf) . (149) 

Comparing fll49p with f l40cp and taking into account that Ps ~ (9 {fib — fi^b) we conclude 
that Cs ~ C ((/i& - fJ'ob)'^). 



5.2.2 Vector sector 

The equations of motion for the vector modes read 
1 d f - 1 dSg^ , r -2 , 



rb'^ dr \ r dr 

+ rCZ'-lf (('-'''' - ''''''' - 3) Go + 2r (-1 + r'b') ^) + Wa^N, {l + 2^ 

(150) 

The boundary conditions we impose on 5(7^ are that they are finite at the horizon and that 
A^'* doesn't get corrected at the boundary, i.e. \\mr^oo9^i = OO Note that using fl37|) and 
fIgU]) we find that 

Tib 1 

As was the case for the scalar modes, it is difficult to solve (11501) explicitly. In what 
follows we will study the behavior of dg^j, close to the phase transition and argue that the 
diffusion coefficient /t defined in fl5^ is continuous across /iq. Close to the phase transition 
we can use (I144p and (11450 together with 

g = l + 0{fib- fiob) (152) 

to solve (11501) for Sg^. We find that the leading order solution to Sg^ is 

5n = J_ (vr - 2 arctan(r6) + 2 ln(l + rb) - ln(l + r^b^)) P^^d^^ 

-u'^d^u^ + O {{fib - fi^b)) . (153) 



r2fe2 + 1 



^^Had we been interested in second order transport coefficients we would have needed to set Unirnex) g^i 7^ 0. 
See [4] for details. 
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Using 



4{J,) = lim {G, + 6G,) + (-3"*) , (154) 

r— >-oo 



we find that T^, the viscous correction to J^, receives a contribution of the form 



= -^P,'9,f + Oifib~ f^ob) , (155) 



from which we conclude that 



7tT 

K=—^ + 0{fib- fiob) . (156) 
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The leading term for k agrees with the value obtained in the normal phase when the charge 
of the black hole is taken to be small [Ml US] . Thus, the diffusion coefficient k is continuous 
across the phase transition. 



5.3 Viscous superfluids at order q 

Our analysis of superfluids in section 12.31 revealed that within our approximation, the stress 
tensor has two possible shear viscosities: the familiar shear viscosity, rj, associated with a 
d{^iUu) term and a shear viscosity, 77^, associated with the gradients of the relative superfiuid 
velocity di^^n^). Standard arguments due to Landau [22] which we presented in 12.31 show that 
?7s = 0. These arguments rely on an assumption of the form of the entropy current and have 
been known to fail in certain subtle cases [T^ [TB| [TB] . The goal of this section is to compute 
?7s and 77 holographically. We will show that 77^ = as expected, and that the shear viscosity 
of the normal component satisfies the universal relation (JT)) as was already discussed in [65] . 

As in section 15.11 we begin by decomposing the order corrections to the metric into 
scalars 6s^'^\ 6f^'^\ a vector ^v/^-* satisfying u^dvli^^ = and a traceless symmetric tensor 
SirliJ which is orthogonal to the normal velocity u'^SirlJ = 0: 

- r^df^^'^u^dx^'u^dx" + 2r^5Vll\^)dx^dx'' + 5nf} dx^ dx" . (157) 
Since the viscous corrections to Tfj_u can come about only through the tensor modes (57r/ti;, it 
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is sufficient to compute T2z- Using ( |7T]) . we find 

kIt,, = -^^23 + 2 lim /5vrg) + ^ Jim r^' {f^'^ - S^s/^s^'^) '^23 + ^ }imna^,+ (l^^) 

(158) 

Here 5a 3/2 is the Kronecker delta function. It appears because limr_^oo r^s'-^^ is non vanishing 
only for A = 3/2. 

(2) 

Since the equation of motion for ^vr/ji/ decouples from the other components of the metric 
we will focus exclusively on it. It is convenient to define 

57r(2j = (57r.(r)a;, + (57r.(r)a^, (159) 

where cr^^ and o"^^ were given in (!29|) . Using ( I159P the equation of motion for SirliJ can be 
decomposed in twoF^ 

A(.(,V-1))1..,)^-.*1m (160) 

i - 1)^*"") = -i - (1 - Ty') ■ 

Integrating fll6ip once and requiring that SiTn is well defined at the horizon, we find that 
lim r'STT^ = -] lim {f^^^ - S^y.s^'^) + ^^^-^ + • (162) 

r— >co 4 r— )>oo ^ io V terms / 

Inserting fll62p into 01580 and comparing to fl33|] we conclude that 

2 1 3/(2) (1/6) ^ ^ 

The (9(g^2) corrections to the shear viscosity given in (11630 are expected. Once f^"^^ 7^ the 
location of the horizon is shifted to 



1 _ f('\l/ b) 
b Ahq^ 



Inserting (I164p into the Bekenstein-Hawking formula for the entropy density, one finds that 
the corrections to the shear viscosity in (I163P exactly cancel the g"^ corrections to the 
entropy density coming from the shift in the location of the horizon (I164p so that the ratio 



^^Equation (|16ip can be read off the q ^ term in an appropriate large q series expansion of equation (27) 
in EQl. 
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of the shear viscosity to entropy density retains its universal value ([T]) up to order 0{q^^). 

To obtain 77^ we integrate (11601) once. Requiring that ^vr^ is finite at r6 = 1 together with 
(11061) implies 

lim r^Sns = — lim rj . (165) 
Inserting (I165p into (11581) and comparing to (l33l) we obtain 

Vs = 0. (166) 



5.4 An explicit solution 

We now turn to solving (I135p and (I136P near the phase transition for the special case = —4 
discussed in [JH] and reviewed in section 14.51 As in section 14.51 we formally expand 6G and 
Sp around /i = /iq = 2/6, 



5G = ^5GE,'"^(r6)(/i6-2)" 

r __ (167) 
5p = ^5p(2-i)(r6)(/i6-2)(2"-i)/2 



n=0 



We will soon see that 5p^~^^ 7^ which implies that 5p diverges close to the phase transition. 
Thus, our analysis is valid only for when the gradient corrections to the various fields are 
smaller than the relative magnitude of the chemical potential, df^u'^ <^ {pb — 2). Since the 
equations of motion for the delta'd quantities (I135P and (I136P are linearized versions of their 
non delta'd counterparts (^^, inserting the expansion (I167p into ( 184|) yields a set of equations 
very similar to those obtained for Gq and p which were described in section The solution 
for the leading order terms in (I167p is given by 



^yU + ^J (168) 



SG^o)(x) = — d - ^^^^ ^ ^^^^ —d N^h 

3(l+x2)'''''' 49x2(l + x2) 

We also need the solutions for 8p'^^ and 6G^q \ but they are too lengthy to repeat here. 
To solve the equations of motion for the vector modes we insert the expansion 

00 

5g, = Y.^gf\rh){ph-2f (169) 

n=0 
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into fllSOp and collect terms with equal powers of fib — 2. The first couple of regular solutions 
which vanish at the asymptotically AdS boundary are 

59^n\x) = ^ - 2arctanx + 21n(l + x) - ln(l + x')) - 

^ All i 1 + x"^ 

5gf\x) = -^^^^YT^ " ^ ^^^^^^^ + 2 1^(1 + ^) - 1^(1 + - 2) ^/^^ ~ a+^T """^""" 

(170) 

The expression for Sg^j^^ has the expected form discussed in section 15.2.21 which ensures that 
K, is continuous across the phase transition. 

With the solutions 6p^~^\ 6G^^\ Sp^'^\ and 6G'^^ in hand we can compute the correction 
to the Josephson condition 

u^d,<P =-p- Jim = -/^ + + ^^23 - 3004 In 2 ^^^ _ ^^^^ ^^^^^ ^^^^^ 

Comparing this to the third line of (133|) and using flllSp we obtain ([3]). As expected, (3 has 
a power law behavior close to the phase transition. Using 

4\{0^)\ = Ca lim (p + 6p) + 0{d') (172) 

r— >oo 

we find that the viscous corrections to call them 6\{0^)\, take the form 

2.,/^ M dV^CAduN" 70993 - 116496 ln2 , , ^, . , , ^.2x\ 

Finally, using f ll70p together with fll54p and ( 15^ we compute 



— {ph-2) + 0{{ph-2f) (174) 



which gives us ([2]). 

The leading order expressions for n and C3 can be obtained from the Kubo formulas (jH] 
and ( l45l) and the analysis of |48] . We present this analysis in appendix |Bl 
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6 Discussion 



In this work we have constructed a slowly varying, spacetime dependent, holographic su- 
perfluid and computed its transport coefficients. The resulting stress-energy tensor and 
charged current of the gauge-theory superfluid ffis nicely with the relativistic version of the 
Tisza-Landau two fluid model. Apart from the shear viscosity and diffusion coefficient, the 
conformally invariant two-fluid model allows for another transport coefficient which we de- 
noted by C3. In section 12.51 we saw that Cs? together with k, and rj, is responsible for the 
attenuation of second sound. Another physical role of C3 can be understood by considering 
entropy production. With the assumptions specified in 12.31 the divergence of the entropy 
current can be read off of (138|) . It is given by 

d,J^ = I (d, (p^nnf + nP^"" {d,^) {d.^) + f^.^^n^ . (175) 

Thus, for n'^ ^ 1, is responsible for entropy production due to compressibility of the 
superfluid component. Indeed, expression (11751) can be compared with the divergence of 
the entropy current in a normal (non superfluid) relativistic fluid which is not conformally 
invariant, c.f. |22]. 



d,J^' = I {d,un' + ^P'-" {d,^) (5.^) + ^a;^d^u'' (176) 

where (i, the bulk viscosity, plays a role very similar to (3. 

The power-law divergence of (3 near the phase transition that we observed in our holo- 
graphic computation does not imply that entropy production diverges there. While ~ 
O {{fib — /io&)~^), the superfluid density vanishes, Ps ~ (9 {fib — Unb). The extra power of Ps 
in (11751) guarantees that dfj_J^ vanishes close to the phase transition!^ A similar observation 
can be made regarding the dependence of on nf. The boundary value of is independent 
of Newton's constant. Since (3 appears together with then ps ~ (which follows from 
(IIMD) implies that Cs ~ ^^s- 

We note that too close to the phase transition one expects the hydrodynamic approx- 
imation to break down since the Landau- Ginzburg potential becomes approximatety flat. 
Our expansion near Tq is then valid only when the hydrodynamic fluctuations are smaller 
then the relative chemical potential. This fact also came into play when we discussed the 
expansion of the bulk scalar near Tq in section 15.41 



^®Wc thank D. Son for pointing this out to us. 



44 



The analysis carried out in this work is vahd for configurations where the charge of the 
scalar field is large. From the structure of the perturbative expansion, it seems likely that 
will diverge near /io for smaller values of q. As long as the boundary value of 6G^ does not 
vanish close to the phase transition, the fact that vanishes there implies through fl40cl) 
that (^3 must diverge. We have no argument for having 77^ = in the fully backreacting case. 
We feel that the result fll66p valid at order together with fl55|) make it unlikely that 

rjs will be q dependent. 
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A Bulk to boundary mapping 

The prescription for mapping bulk fields to boundary expectation values can be found in 
sections 14.11 to 14. 3[ In what follows we collect these results for ease of access. The energy 
momentum tensor of the boundary theory is given by 

i4{Tt,u) = - hm ( - K-f^^ + 87^^ - ^Tmi^'^aIV'P 

- ^IfiuffiA {ip*n"'daip - {il)*n^dyi) + ip*n^df,ip) + c.c.) ^ . ([71]) 

The norm of the scalar follows from (173|) . 

4 1 (O^) I = Ca hm [p + 6p) + 0{d^) . 
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The expectation value of the gradient of the Goldstone boson is given by 



= _ hm (Go + 5G) . (dH]) 

r— >-oo 

The expectation value of the current follows from (176|) 

kI{J,) = Jim {G^ + 6G,) + [^^^ . (llMj) 
Using our definition of the charge density in section [2] we find that 



KgPt = - hm r^Go + ( * 



vci'gcnt 
terms 



(unD 



9 / possible \ 

- lim r g + [ divergent 1 . 
r^OO V terms / 

B Two Point Functions from Gravity 

In this appendix we work out the leading term for n and ^3 using the Kubo relations in 
section 12.41 Most of the Greens function needed for carrying out this computation have 
already been computed in |38]. We briefly go over the notation of |38] and compute the 
remaining Greens function. 

We start with the action fl64p and work in a weak gravity (or probe) limit g ^ 1 in which 
gravity decouples from the scalar and gauge fields (or Abelian-Higgs sector) [S]. In this limit, 
a solution to Einstein's equations is a black brane: 

ds' = ^ (-f{u)dt' + dx' + dy' + dz' + ^) , (177) 

where f{u) = 1 — {u/uhY- At m = > 0, there is a horizon, while in the limit m — > 0, the 
space asymptotically becomes anti-de Sitter. Since we will use linear response to compute 
the transport coefficients it is sufficient to work in the Fefferman-Graham coordinate system. 
In (!79|) we used an ingoing Edington-Finckelstein coordinate system. For convenience, we 
will set = L = M/i = g = lin what follows. We also work in a gauge where /I5 = and 
consistently set ip to be real. 

When m? = —4 the near boundary expansion for tp and takes the form 

iP = {^^''^u' ln{u/X) - {O^y + ...) , (178) 
^M = 4'^ + ^(^>' + -- - , (179) 
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where we have introduced a UV cutoff A ^ 1. In this coordinate system, in place of f ll72p 
and f ll54p we have 

(J^) = hm , — -- , (O^) = hm — — , (180) 



where gliJ = hm^j^o u g^i, is the Minkowski metric. Since we are working in a gauge where ip 
is 



is real, we can identify Af^ with the chemical potential fi and Af^ with a superfluid velocity 
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As described in section we can solve for ip and A^ near the phase transition fi/T = 2tt. 
Following the notation of [H] instead of using ^/T — 2tt as a small parameter, we use 

e = -V2{0^) . (181) 

At equilibrium in the absence of a superfluid velocity, we have the near boundary expansions 

^/2ijj =eu^ + 0{u'^) (182) 



At = {2 + 5fi2e^ + 5/14 6^ + 0(6^)) - (^2 + Q + 5/i2^ 6^+ 
-5 + 6 log 2 ^ \ ^4 ^ \ ^ ^ 



;i83) 



1152 



where 



1 . 253 71og2 , 

(5/i2 = — , (5/i4 = 1 — . (184) 

^ 48 V 55,296 1152 ' ^ ' 



To compute the Greens functions we consider the fluctuations of the fields: 

SAt{u,t,x) =at(M)e-'"*+'^" 

M^(M,t,x) =a^(M)e-^^*+^^^ (185) 
Si;{u,t,x) =4j{u)e-''^'+'''yV2 . 

The Greens functions can be computed from the coefficient of the term in a near boundary 



^^There is a minus sign discrepancy between this identification and the one in [6]. Here we define ^4^(0) — 
whereas in [6], was defined to be the phase gradient of the scalar. 
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expansion of the fluctuations (see [H] for more details). We find: 

- ^ [iSik^ + 7ie^ + USe^u - miu^ + 8P{5ie^ + 12a;)] 

^ [I2k^ + (Te^ - 120iuj)uj + 3k\e^ + 16(1 - i)uj)] 

+ ^ [A8k^ + + 12(2 - 3i)e'^u - 96(3 + Ai)u^+ 

+16A;2(e2 + 6(l-2i)w] (186) 

4A;2(7e2 + 72(1 - i)w 











uk 








G*« 









^oo ^ ^oo =^(4A;^ + - 32za;) 

where the pole structure is given by 

V = 960a;^ + 56z(12A;2 + e^)^^ - 12(16P + 3e'^)k'^u - z(48A;^ + Wk'^e'^ + e^)P . (187) 

Expressions fll86p and f ll87p are the leading order results in a limit where ~ /c^ ~ ^ 1 
are all small and of the same order. The complex conjugate Green's functions involving 
can be obtained by sending u — )■ —u, k — —k, and i — —i. These Green's functions obey 
the Ward identities -wG" + A;G*^' = 0, -wG*^' + kG""^ = 0, -wG^* + kG'^^ = (O^), and 
—coG'^^ + kG'-^^ = — (O^). From the two-point functions involving O and O, we can construct 
a two-point function for the phase of the condensate: 

Q<p^ ^1 ( ^oo ^g'^- G^^ - G^^ 



2e 



4i 

(48fc^ + + 16A;2(e2 - Qicj) - 32ie2a; - 96^^) . 



e2p 

By direct computation, we find the following Kubo relations 

lim hm — Re G" = lim lim Re G^^ = , (189a) 

u)^Ok-fO k^ w^Ofe^O 4 

lim lim ^ Im G" = lim lim ^ Im G''* =1 , (189b) 

oj-i-O fc-i^O k'^ uj^O k^O k 

lim hm Re G'^'^ = - — , (189c) 
w->ofc^o 14 

52 1 

lim lim u Im G'^^ = . (189d) 



48 



Comparing with (H6|) , (H5|) , (HHj) , and (H9|) , the right hand side of these Kubo relations should 
correspond to — Ps//^) i^/T, —djji/dpt, and (^3 respectively. Restoring factors of K5 and T, we 
find that they do agree at leading order with (11221) . ([2]) and ([3]). Since f ll89bl) is independent 
of the value of ps, the relation implies that k is continuous through the phase transition. 

Another way of extracting k and is through attenuation of fourth sound. The solutions 
of P = give poles of the Green's functions from which we can infer the dispersion relations 

a; = ± — = h..., uj = , (190) 

2a/T4 196 2 ' ^ ' 



for T < Tq and T > Tq respectively. These expressions should be compared with (!62|) and 
( 163|) . Note that Ps = and dpt/dp = 2 for T > Tq while for T < Tq, these quantities are 
given by (I189ap and (I189cp . If we read off the value of k from the T > Tq result, we find 
agreement with (]189bp . If we then assume that n is continuous through the phase transition, 



we can deduce the value of Cs and find agreement with fll89dp 



C Formulae for the Explicit Solution near Tq 

This appendix provides the first few terms, in an expansion in {p — po)b, of the static 
background solution near Tq described in section 14. 5[ For the scalar quantities Gq and p we 
have 



x^(x^ + 1) 

(4). _ 309x^ - 56a;2 - 253 , (72^2 -96) In 2 , 24:xHn{l + x- 



4a;2(a;2 + l)2 ^ x^{x^ + 1) x^ + 1 (191) 



1 + x^ 



, 253x2 + 157 28y31n2 4^3 In (l + 4) 

0^ '(x) = 1 

^ ^ ^ 4v^(x2 + 1)2 x2 + 1 a:2 + 1 
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For g, the first few terms are 



g(^\x) =1 



a;2 + 1 



9^'\x) =6 (hi, {-x')+U2 {l-x')+U, (^)) 

289x^ + 229 WM ^ 2\ 12x2ln(l + 34in(2) 



- 61n(2)ln (l - x') + -r^-^^ ; ^ »^ ^^H-; + 31n^(2) 

3 (2 In {x^ (l - x^)) - In (x^ + l)) In (x^ + l) , 



4(x2 + l)' 



For s^^^ and Z^^-*, the first few terms are 
=0 

s^'~'\x) = ^ 



48x^ + 421x^ + 618x2 + 221 112 ln2 16x^(2 + x^) ln(l + x" 



2^ 



3(1+X2)4 (1+X2)2 (1 + X2)2 

4 



3x6(1 +x2) 

4) _ 4 (24x1° + eOx^ - 208x6 - 335x^ - 38x2 ^ 



3x6(1 +x2)3 

64(5x2 - 2) In 2 32(x^ + x^ - 2) ln(l + x"^) 



X6(x2 + 1) x2(x2 + 1) 

Finally, the vector mode corrections to the metric are given by 

7® (a;) =0 
7®(x) 



(192) 



f^-'^i^) (193) 
j(2,2)^^) 4 (-5 + 19x2) 



X4(x2 + 1) 

^^^\x) = - ^^(^) + ^Gf^ - f- \^ ^ (413 + 538x2 + 269x^) + ^60 In 2 . 
X* 2x^ 8(x^ + l)^x* X* 

(194) 

Note that the leading order corrections to the backreacted metric /(2'0)and s*^2>o) together 
with the leading order expression (179!) generate an exact solution to the Einstein equations — 
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the Reissner-Nordstrom black hole. This solution describes the fully backreacted geometry 
at high temperatures when is not condensed. 
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